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The strong coupling effect in the (2+i)-dimensional quantum electrodynamics (QED3) description 
of the 5 = 1/2 Ifeisenberg antiferromagnet is studied in terms of a canonical transformation which 
has been used in the small polaron theory in electron-phonon systems. In the Neel ordered phase, 
we show that the Dirac fermions are localized and the spectral function is of the Gaussian form due 
to the coupling to the longitudinal gauge field. The width of the broad line shape is ~ 3 J with J the 
superexchange interaction energy and the localization length is ~ 2a, with a the lattice constant. 

PACS numbers: 



I. INTRODUCTION 

One of the most fundamental question in the physics of high-Tc cuprates is about the description of the doped 
holes. The simplest situation would be realized in the single hole doped system. Experimentally, angle resolved 
photoemission spectroscopy (ARPES) in the undoped compounds provides us valuable information about the nature 
of the excitations. The excitation spectrum was obtained by Wells et almt^ It was found that the band width is 
^ 2.2J ~ 270meV where J is the superexchange interaction between the copper site spins. This band width is 
much smaller than the tight-binding value of 8t ~ 2.8eV which is determined by the transfer integral parameter, 
t. Furthermore, the observed spectra exhibit a quite broad feature. In Ref^, it is shown that the spectra do not 
have a conventional Lorentzian form. The spectra rather have a Gaussian form. These observations suggest that 
the quasiparticle excitations observed in the undoped cuprates are quite different from conventional Fermi liquid 
quasiparticles. 

In order to understand ARPES experiments in the undoped compounds, first we need to figure out the elementary 
excitation whose excitation spectrum obeys the dispersion observed experimentally. The next step is to consider a 
mechanism of spectral shape broadening that would arise from coupling to a boson mode. Concerning the first point, 
it was shown that the slave-fermion approach to the t-t'-t"-J model reproduces the dispersion>^ (It is known that 
the simple t-J model fails to reproduce the dispersion along (7r,0) to (0, tt).^'^) However, it has been shown-^ that 
the coupling to antiferromagnetic spin-wave modes does not lead to a broad line shape. Recently Mischenko and 
Nagaosa^ considered a coupling to an optical phonon. They numerically summed over Feynman diagrams including 
vertex corrections for phonons. It was argued that the quasiparticle in the t-J model is in the strong coupling regime 
that leads to a broad line shape. 

In this paper, we consider another possibility. We consider the quasiparticles in the tt-Aux phaseiSi and the effect of 
coupling to U(l) gauge field fiuctuations. A similar problem within the slave-fermion approach to the t-J model was 
discussed by Auerbach and Larson^ It was argued that the single hole forms a small polaron due to the coupling to the 
gauge field. Here we consider a similar problem based on the 7r-flux phase with dynamically induced mass. The 7r-flux 
phase was proposed by Affleck and Marston-'^^ from a mean field theory of the S* = 1/2 antiferromagnetic Heisenberg 
model based on a fermionic representation of the spin S = 1/2. Fluctuations about the state can be described by a 
U(l) gauge field. The dispersion of the quasiparticle in the 7r-flux phase is in good agreement with the experimentally 
obtained dispersion as argued by LaughliniiS, The effective theory is described by quantum electrodynamics in (2+1)- 
dimension, which is called QEDaiiSii^ 

Although the 7r-flux phase is different from the Neel ordering phase, believed to describe the Neel ordering stat o^'^i^^ 
is the state that includes the dynamically induced mass term due to the coupling to the gauge field. The excitation 
spectra near (±7r/2, ±7r/2) are given by a massive Dirac fermion. Experimentally the Dirac fermion mass is estimated 
to be ^ 1.3J^^ by fitting the experimentally obtained dispersion with the quasiparticle dispersion in the tt-Aux phase 
with the mass term. 

Meanwhile Franz and Tesanovic proposed^^-^ a QED3 theory starting from the d-wave superconducting statei^ The 
dynamical mass generation, or the spontaneous chiral symmetry breaking, is also associated with antiferromagnetic 
orderjiS Aitchison and MavromatosiS argued that the gauge-fermion interaction in QED3 leads to non-Fermi liquid 
behavior. Rantner and WenSfl^Sii discussed changes of line shapes between the normal state and the superconducting 
state based on a flux phase obtained in the SU(2) slave-boson approach to the t-J model^^ with QED3 for the 
description of the spinous. 
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In the Neel ordered phase, there is strong logarithmic potential between the Dirac fermions. In fact, there is no 
low-lying spin 1/2 excitations22i24 because of the confinement associated with this strong interaction.^^ Therefore, in 
the Neel phase we need to perform a non-perturbative analysis. 

In this paper, we apply to the QED3 a canonical transformation, which has been extensively studied in the electron- 
phonon systems, for the effect of the longitudinal gauge field interaction on the massive Dirac fermion. We show that 
the coupling to the longitudinal gauge field leads to the localization of the Dirac fermions and the spectral function 
is approximately given by a Gaussian form. The main assumption here is that a doped hole excite a quasiparticle in 
the TT-flux state. One way to interpret this is to consider the slave-boson approach to the t-J model. Another way 
is proposed in Refill where a half-skyrmion spin texture is created by doped holes. The excitation spectrum of the 
half-skyrmion spin texture is the same as that of the 7r-flux phase. The mass term of the quasiparticles in the tt-Aux 
state is the excitation gap of the half-skyrmion spin texture. The effective action of the half-skyrmion spin textures 
is shown to be described by QED3 by applying a duality mapping. 

The rest of the paper is organized as follows. In SecEl we describe the QED3 action and the application of the 
canonical transformation. In Sec lIIII we calculate the retarded Green's function of the massive Dirac fermion. Section 
IIVI is devoted to the conclusion. 

II. MODEL 

In the continuum, the QED3 action for the = 1/2 quantum Hcisenberg antiferromagne1ii^*i^ is written as 
S = J d?x^ {x) [i^D^, - mCwCTs) i^{x)-^ j d^xF^^F^''' 

= So + Sint + Sa, (1) 

where Z?^ = 9^ + iA^^ is the covariant derivative. (The derivation of the action in the tt-Aux phase is given in 
appendix^) The Dirac fermion fields describe the quasiparticles in the spin system. In the absence of the doped 
holes, the ground state is that all the negative energy states of the Dirac fermions are occupied. While the positive 
energy states are empty. Upon hole doping, hole states are created in the band of the negative energy states. In 
writing down the action, we have assumed that there is the Neel ordering. In the presence of the Neel ordering, 
the mass m is nonzeroii^i^ The dynamics of the gauge field arises from the fermion polarization.^'' Due to the mass 
term the gauge field action takes the form of the Maxwell action. The gauge charge ca is related to m through 
e\ = 3nmcsw '^'^ 

For static particles, the scalar potential (j) = Aq/ca is calculated by the following equation: 

V^^iv) = -e^p(r), (2) 
where p(r) represents the charge density. The potential between a particle and hole is given by 

Vir) = -^\n-. (3) 

Here tq is the short distance cutoff scale which is on the order of the lattice constant. Note that the potential V{r) 
has a form of a logarithmically confining potential. 

We consider the effect of the scalar potential on a dynamical hole. In solid state physics, the Coulomb gauge is 
often used for the analysis of the electromagnetic gauge field. In the Coulomb gauge, the interaction mediated by 
the scalar potential is instantaneous in time. While in the Lorentz gauge, the retardation effect in the interaction 
is explicit. For the electromagnetic gauge field, the speed of light c is much faster than vectors of excitations in the 
sample. Therefore, the retardation effect is negligible. However, for the QED3 theory of the speed of the gauge 
field propagation is Csw that is much smaller than c. The retardation effect can play an important role. To include 
this effect, we take the Lorentz gauge. In the Lorentz gauge, the action of the gauge field is 



/j=0,l,2 

The vector potential obeys the following equation, 

1 



SW 



(4) 



dfA-dlA-d'^A^-l (5) 
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Note that the equation for the longitudinal component of the vector potential is derived from the equation for the 
scalar potential in the Lorentz gauge. Therefore, the independent degrees of freedom are the transverse component 
of the vector potential and the scalar potential. The Hamiltonian for the gauge field is 

where uiq — CswQ in the continuum and we have introduced creation and annihilation operators through Aq^ — 
V e^Csw/(2wg) (^alg,^ + o-qfi^ ■ Here /i = t denotes the transverse component of the vector potential. 
On the square lattice, the action So + Sint has the following formSS 

So + Shit = J dta^J2 {^3 (*) n"dt^j (t) - (t)csw7°^o(i, t)^j (t) - mclJp^ (t) ag^j (t) 

3 

Za I J 

(t) [e-^^(^'+'^^=/^'*V,+ae. (t) ^^(^--^/2^*)^,_<,g, {t)] } . (7) 

The above action contains the gauge field interactions arising from both the transverse part and the longitudinal part 
of the gauge field. The transverse component of the vector potential is important for the calculation of the quasiparticle 
lifetime. Because the scattering process associated with the transverse component is due to emissions or absorptions 
of the propagating gauge fields, or "photons." As demonstrated in Refsi22i2£, the coupling to the transverse vector 
potential leads to a T-linear resistivity law. However, the interaction of the holes with the transverse vector potential 
does not creat a trapping potential, and is described by a self-energy effect. Therefore, the spectrum has the form 
of Lorentzian which is inconsistent with the ARPES experiments in the undoped compounds. Here we focus on the 
self-trapping phenomenon associated with the scalar potential, and neglect the transverse component of the vector 
potential: 

Hhit = Cswa'^ -0,-7°Ao(j)'0j- 

3 

= ^EE^lV'.e-^^A,. (8) 



The interaction i^int contains the terms that contribute to the infrared (IR) behavior of QEDg^^-^^-'^^ and the other 
terms. As is well-known in QED3, the dynamical mass generation is associated with the IR behavior For the 
analysis, we need a non-perturbative scheme, such as a 1/7V expansion with N the number of Dirac fermion species. 
For discussion of the dynamical mass generation, the Schwinger-Dyson equation has been studied. '^^ In the analysis of 
the Schwinger-Dyson equation, it is necessary to figure out the form of the vertex part. For the IR region, one needs 
to properly choose the form of the vertex part for the discussion of the dynamical mass generation. While for the 
other region, the bare vertex form is known to be sufficient for the analysis of the system^ Here we are interested 
in the interaction with the momentum which is not in the IR regime. The IR behavior of the system is understood 
to be associated with the mass generation. Therefore, we introduce an IR cut off qo, which is much smaller than the 
ultraviolet cut off A ~ 7r/a, in the following calculations. 

We consider a canonical transformation:^^ 



H = exp{s)He^p{-s), (9) 



where 



with = a^Q and = a^o, and 



^ = E ^ - ' (10) 




,2 „3 
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The bare vertex is used here as discussed above. Hereafter the summation in the momentum space concerning the 
gauge fields is to be performed for go < IqI < A. 
Using the following result, 

we find that the fermion field tpj is transformed into 

e'lPje^' = XjiPj, (13) 

where 



Xj = exp 
The transformed Hamiltonian reads 



2 / \ 3/2 



X! 7~ exp(iq-Rj) (aLg 



(14) 



H = Ho + i/int, (15) 

where we have omitted the band narrowing factor X'^-Xj±ae^^, in the hopping term. (The estimation is similar to that 
of eas. H38|l and H39|) below. See Appendix IbI) Hmt describes the logarithmic interaction between the Dirac fermions. 
In the following, we consider the effect of the longitudinal gauge field on the Dirac fermion spectrum given by Hq. 

III. CALCULATION OF THE RETARDED GREEN'S FUNCTION 

The retarded Green's function may be given by 

Gr (z, j; t) = -id it) (^A, (0) (0) V. it)) (18) 

We separate this quantity into two terms, 

GRii,,r,t)^Giii,j;t) + G2it,j;t), (19) 

where 

Giii,,r,t) = -tdit){Ait)i^jiO)), (20) 

G2(i,j;t)--z0(i)(^,(O)V'.(i))- (21) 

In the following calculations, we only consider the Dirac fermions with the spin up component and residing near 
(7r/2,7r/2) in the momentum space. We redefine the Dirac fermion field tpi^) so that it represents the spin up Dirac 
fermion and describes the components near (7r/2, 7r/2). For the 7 matrices, we may use the reduced form of 7 matrices: 
7^^ — > T3, 7^ — > iT2, and 7^ iri with the Pauli matrices because we only consider the two components of the 
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Dirac fermion fields. (The calculations below are equally applied to the Dirac fermions with the spin down component 
and/or residing near (— 7r/2, 7r/2).) For t > 0, by performing the canonical transformation we have 

iGi{i,j;t) = (ViWV^,(0)) 

Tr e-'^"»e^'>*XiiPie-'"°*^jXj 
Tre-0"o 
eA 



exp 



fe.fe' 



g UJg 



-270 (i sin kiji+i sin k2J2+m)t^^^^^ ^ik-Rt ^—ik'-R 



X exp 



(22) 



where we have used 



e'^°Vfce"*^°* = exp [i (-n sin fci + T2 sin A;2 + rsm) Vfe- (23) 
(Hereafter we set a = 1 and Csw = 1-) The matrix in the exponential is diagonalized by the following unitary matrix: 

1 / Ek +m 



Uk 



sin ki -\- i sin k2 
^2ek (£fe + mj y-^'i-'o-k-i + i sin A;2 £fc + m 



£;j + m — sin fci — i sin fc2 



V2£fe {ek + m) V ^i'^ " * ^™ ^2 £fe + m 
with Cfe = V sin^ fci + sin^ + m^, and then, 

iGi{i,j;t) = (exp 



eA 



xexp 
1 



2f2 q UJq 



X 



fe 



[l-/(£;fc)] e-''"* 







jji^ik-iRi-Rj 



Using the following formula, 



,Aa^+Ba Ca^+Da _ ^{AB+CD)+BC JA+C)J {B+D)a 



(24) 
(25) 



(26) 



(27) 



we obtain 



exp 



= exp 



= exp 



eA 



2^ q l^q 



exp 



2^ q 



(28) 



£A "S^ J_ /"l _ pioJ,t -iq-{Ri-R 

2f2 ^ u;3 V' 

^ ^ sin [wgi -q-{Ri- Rj)] 



exp 



g2 

"jf E j;;^ (1 " ['^9* - g • (i?i - i?j)]) '^K) 



20 w„ 

q 9 



X exp 



ll E i coth ^ (1 - COS [a;,i - 9 • (i?. - i?,)]) 



(29) 
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Since ek > m — J , and we are interested in the temperature region where PJ ^ 1, we may assume that /(efe) ^ 1. 
Thus, 



iGi t) ~ exp 

X exp 
1 



— sin [LOqt - q ■ {Ri ~ Rj 



2rj ^ w„ 

q 1 



ik-(Ri-Rj) 



1 



Performing similar calculations, we obtain 



- sin — i sin ^9 



2£fc \ ~ sin fci + i sin fc2 et — m 



iG2 {i,j;t) ~ exp 



217 



X exp 



^T.^ COth ^ (1 - COS [coqt - q ■ (i?,: - R,)]) 
j_ / 



Finally, we obtain 



iGR{i,j;t) ~ exp 
1 

X — 



fi ^ i. COth ^ (1 - cos Kt - g • (i?. - i?,)]) 



q 1 



E 



^-iekt^ik-{R,-R,) p (ij-t) 



with 



f{i,i\t) = exp 



^ sin (wgt) cos [(7 • — Rj 



2n ^ uj„ 

q 9 



(30) 



(31) 



(32) 



(33) 



where we have assumed that [m + efc)/2efc ~ 1, (m — ek)/'2^k — 0, and (sinfci ± isinfc2)/2ei; ~ 0. Because for the 
small hole concentration, k/m <^ 1. Note that the factor. 



(34) 



denotes the wave function of the free propagating particle. 

The second exponential factor (= /) in the last line of Ea. (|29|l leads to a localization effect of the Dirac fermions. 
/ is approximated as follows. 



f dq f de-^COth^il- COS [uJqt-q-{Ri-Rj)]). 

Jq Jq ^q 2 



87r2 



(35) 



From numerical estimations, we see that / has the Gaussian form as shown in Figs. and |21 We find that 

/~exp(-1.8125i2- 0.90625i?2). (36) 



To reproduce this Gaussian form, we make a following approximation: 



I ~ 



d^q 1 (3ujq 



+ {q ■ [R, - Rj)Y 



{Ri - Rj) 



(37) 



7 



0.8 



0.6 



0.4 



0.2 



t=0.0 
t=0.6 
t=1.0 




where 

-12 



A^ = 4/^^coth^ (38) 
J {2ny uj„ 2 



At low temperatures, A is 



1 



^ q (iu 
dq — coth 

qo ^1 2 



0? q q ^, P^q 

— —^^coth— ^ 
27r)2 uj^ 2 



(39) 



(A - go) 



n=l 



-nl3qo 



-n/3A 



n/3 



Using the values of Zc = 1.17, which was estimated from quantum Monte Cairo simulations,^^ and m = 1.3J, which 
was estimated experimentallyfi^ we obtain A ~ tt/o = \/2ttZcJ — 5 J, e\ = Sttto ~ 12 J. Thus, 



A ~ 3 J. 



(40) 
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This value is in good agreement with the value A ~ 2.7J, which is obtained from Ea. (|36|l . Therefore, the above 
Gaussian approximation is justified. Similarly H. is evaluated as 

^:.^y^^d,cothY-^A, (41) 



£~ 2a. 



(42) 



From Ea. l|36|) . we find I c:^ a. These result suggest that the particle is localized with the localization length oil = 0{a). 
The function f{i,j',t) is associated with the wave function rcnormalization. The argument of the exponential is 



1 



A 



90 



-Jt\n — < In — 



A 



90 



(43) 



(The value of t that gives the saddle point in the frequency integration may be used in the calculation of Gii{k,u)) 
below.) As a rough approximation, we take f{i,j',t) as a constant C. Thus, 



iGR {hj-t) = exp 



4 



\Ri — Ri 




(44) 



Performing Fourier transformation, we obtain 

_ 87r5/2 ^Cexp[-(c^-£fe)VA^] 




iGR{k,uj) 
The spectral function is 

-llmG.(fc,.).Z^f-P[-('^^^^)VA^] 





C* exp[-(a. + efe)VA2 



exp [-(t^ + efc)VA2] 



(45) 



(46) 



where Z is the wave function rcnormalization factor given hy Z — 47rReC. Thus, the spectral function has a Gaussian 
form. This is consistent with the experiments.- The value of the width A ^ 3 J appears reasonable compared with the 
experiments. The main drawback in the present calculation is that we cannot discuss anisotropy in the momentum 
space observed in the experiments. Because we have started from the continuum theory. In order to discuss anisotropy, 
we need to formulate the theory on the lattice. 



IV. CONCLUSION 



To conclude, we have studied the effect of the longitudinal gauge field on the quasiparticle spectrum in the tt-Aux 
phase with the dynamically induced mass. A canonical transformation is applied for the analysis of the strong coupling 
effect between the massive Dirac fermions and the longitudinal gauge field. We have shown that the Dirac fermion 
propagator has the strong decay form that leads to the localized state. The spectral function is approximately given 
by a Gaussian form. The width of the broad spectral line shape is ~ 4 J and the localization length is ^ 3a. 

Although we have carried out the analysis in the Neel ordering phase, in the disordered phase we expect a qual- 
itatively different behavior because there is no strong logarithmic potential that leads to confinement. It would be 
interesting to study the delocalization transition associated with the Neel order destruction. 
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APPENDIX A: DERIVATION OF QED3 ACTION IN THE tt-FLUX PHASE 

In this appendix we derive the QED3 action in the tt-Aux phase^" with the dynamicaUy induced mass m. The 
S ^ 1/2 antiferromagnetic Heisenberg model is given by 

i7 = J^S, -Sj, (Al) 

(ij) 

where the summation is taken over the nearest neighbor sites. The spin 1/2 can be represented by fermion fields, 



where f ja and /J^ are fermion annihilation and creation operators at site j with spin a, respectively. The fermion 
operators satisfy the constraint, fj^fja = 1- In this fermion representation, the Hamiltonian reads 

^ = E E fLf^^fjpU + * E ( E 4/^" - 1 ) + ^on^t- (A3) 

Here Xj is a Lagrange multiplier to impose the constraint. Following Affleck and Marston^iS we take the mean fields. 



We set xi = Xj.j+e^, X2 = Xj+e^.j+cy, X3 = Xj+e.^+ey.j+ey, and X4 = Xj+ey,3- Now the HamiltoniE 

^ = E E (^l-^l+i,^-^^'^ + ^4/|+y,<x/j'^ + ^-C-) - ^ E {xlf]+i:,J]a+X2f]+y^j3a + h.. 

j^even a j^odd 



+ ^(|XlP + |X2p + |X3p + |X4n+. 



Performing the Fourier transformation 



with X)k' boing the summation over the half of the Brillouin zone, we obtain 



k 



where the two-component spinors fek and fok are /ek = (/k + /k+Q)/\/2 and /ok = (/k - /k+Q)/\/2, and = 
In the Neel ordering phase, there is the staggered magnetization term 

J2 msti-iyf^Tsfj - mst^' (/i^s/ek - flk<J3fow) , (A8) 



where 173 is a Pauli matrix in spin space. Following Affleck and MarstonjiS we take Xi = X3 = Ixl and X2 = X4 = *lxl- 
Then, 

H = y'(A A)( n u Z'""'- .^ ^IxKcosfci+zcosfc^) W /ek U eonst. (A9) 
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The low-energy Hamiltonian near {±tt 12,11 /2) points is 

^ — 5Z ( -^elk /oik /e2k fo'. 



o2k 



T3 
-T; 



mstCTz - J\x\ki 



Tl 
-Tl 



-J\x\k2 









:)] 




(AlO) 



Here /e,oik describes the field near (7r/2,7r/2) and fe,o2k describes the field near (— ■7r/2, 7r/2). The matrices Tj are 
Pauli matrices. Note that in this representation (— 7r/2, — 7r/2) and {tt/2, —tt/2) are equivalent to (7r/2,7r/2) and 
(— 7r/2, 7r/2), respectively. The Lagrangian is given by 



idt - Jlxlidx 



Tl 
-Tl 



J\x\idy 



T2 
T2 



T3 
-T3 



(All) 



where V'U'^) = (l/v^)X]fc' *5^P ("^l' ■'')( /Jik /oik /e2k /o2k ) ■ P^''^°™^^S 



V'(r,t) 



1 

iT2(Tl 



4,{v, t), 



we obtain 



where ip — i/'Wi 



Setting J\x\ 



L = / (friPir,t) [ij^dt - Jlxh^idi ~ J\x\l^id2 - mas] ^(r,t), 



70 



we obtain 







f iT2 












71 = 


[ 


-iT2 j ' 


72 = 





-ITi J 



(A12) 



(A13) 



(A14) 



(A15) 



where mc^^ — TOst- Fluctuations about the mean field state are included by a U(l) gauge field, ^^lifliM This is done 
by replacing with the covariant derivative = + iA^. 

The action of the gauge field is obtained by calculating the fermion polarization function as explicitly shown, 
for instance, in RefflA. Since the fermions are massive, the action has the following Maxwcllian form: 



1 

4^ 



(A16) 



The coefficient 1/e^ is dependent on the mass through = 37rm»i^ 

We derive this action in the imaginary-time formalism. If we take the imaginary time as it — r, and define 
2^3 = CswT, then the action takes the so-called Euclidean action: 



(A17) 



j d^x^p (x) [7^ {d^ - iA^) + M] ^ (x) , 
where M = mcsw The 7-matrices are changed as follows. 



i'-f^ 71, ^7^ 



72, 7 







73- 



Their explicit forms are 



71 



-T2 
T2 



72 = 



-Tl 
Tl 



, 73 



T3 
-T3 



(A18) 



(A19) 



The Dirac fermion Green's function is 



One obtains, 



{iksT3 — ikiT2 — ik2T\ + M)~ 



1 



Thus, the Green's function is given by 



Gk = 



-iks + M ki+ ik2 
M2 + A;2 V -^1 + ik3 + M 



{ik^^j^, - M) . 



k'^ + M'^ 

Now we integrate out the Dirac fermions fields. After Fourier transformation, the action reads, 

S = '^'4'k [{ik^-it, + M) 6k,k' - h^An {k - A;')]V'fe'- 



fe.fe' 



Integrating out the Dirac fermions yields 

S = -Tr In (-G^i) - Tr In [1 + iG^M 
= -rrln(-Gfc^) - Tr 



iG^^A^ - \ {iG^^A^f 



The effective action of the gauge field is 



5(2) 



(27r) 



The polarization function 11^,^ (g) is 



/(fik 



Here tr is taken over Dirac fermion components and spin space. 
In order to calculate 11^,^ (g), we use the Feynman trick: 



By making use of the Feynman trick, we obtain 



— = [ dx 

ab Jo [a (1 - X 



{l-x) + bxY 



dx ■ 



xtr 



(k + xqf + a; (1 - x) g2 + M2 

(«fc^7p - M) 7^ (i {k + q)^ 7m - 7. 



dx ■ 



xtr 



(27r)^Jo [P + a;(l-a;)g2 + M2]^ 
{i {k - xq)^j^ - m) 7^ (z (fc + (1 - x) g)^7^ - m) 7, 



To calculate the terms including 7-matrices, we use the following formula: 



1 



Thus, we obtain 

n,,. (q) - -8 / / dx 



{2TTy Jo [k^ + X {1 - x) + ]VPY 
X [~2k^K + 2x{l~ x) {q^q^ - q^5^,^,) +[k^ +x{l~ x) + M^] 5, 



21 



where we have used trl = 8. 

As a convenient calculation scheme, we use the dimensional regularization»^ 



(27r)^r(i^/2) 7o Jo [k^ + x{l-x)q^ + KPf 

d^k [P + x{l-x)q^ + M2] 6,,^ - 2k^k^ 



{2tt)^ Jo [k^ + x{l-x)q^ + A-Pf 



The integrations in the first term is calculated as follows 

l,D-l 



dxx{\~x) I dk ^ 

Jo [fc2 + x (1 - x) g2 + M2]2 

TT l~D/2 f\^^^^_^^\^M^^x{l-x)q^]'"^'^ 



2 sin {7^0/2) Jo 
Setting IP = 3, we obtain 

dxx{l- x) [M^ + x{l-x)q^^ - 



2q' ./f + UP 

The second term is rewritten as 

d^k f\ \k'^ +x(l-x)q^ +M'^]5.^-2k^k^ 
' dx ■ 



{2TTr Jo [k^ + x{l-x)q^ + A'PY 

[ [\ 1 ^ [ d^k f\ k^k, 

'""J {2nfJo k^+x(l-x)q^ + M^ J ^27rf Jo [P + x (1 - + M^f 

The integral over k in the second term in the right hand side is carried out as follows 

{2nf [k^ +x{l-x)q^ + NPf 

27r-D/2 i,D+l 



-(5^^ / dk ■ 



{2iTyDT{D/2) Jo [k-^ + x{l-x)q^ + APY 
7r^/2+i [x{l-x)q' +AP]'''^'^ 
" 2(27r)^rp/2) sin(^Z?/2) ^^"^ 

The first term in the right hand side is, 

d^k 1 

{2-k)° fc2 + a;(l-x) q^ + A'P 



dk ■ 



{27rfT{D/2) Jo k^ + x {1 - x) q^ + ]VP 

^D/2+l [a; (l-2;)g2 + M2] 
{2nfT{D/2) sin {nD/2) ' 
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From these equations, we find that the two terms in the right hand side of Eq. ljA36p exactly cancel each other. 
As a result, we obtain 



n^i' (<?) = - [Q Sf,„ - qfj^qi, 

TT 



M q^-AM"^ _i q 

TT^ ^ r^5 sm , = 

4g3 y^M^lAP 



Expaiding the right hand side with respect to g, we obtain 



10 \M) 560 \M, 



On the other hand, the Maxwell term has the following form 

Sm = TTT X! W'^i^^ ~ (-9) (9) 



(A39) 



(A40) 



(A41) 



Thus, we obtain = inM . 



APPENDIX B: THE BAND NARROWING FACTOR 



In this appendix, we evaluate the band narrowing factor omitted in ea. p6|l . The canonical transformation (PJ 
introduces the following factor in the hopping term. 



X\X, = exp 



eA 



-q a, I e 



iq-Ri 



exp 



SA 



y- 

^ ,,,3 



a_„ — flo e 



By taking the thermal average with respect to the bosons, we obtain 



A,/A,)=exp 



2Vt 



E ^ (1 ~ cos \q ■ (i?. - R, 



)]) 



The calculation is similar to eas. 1)38(1 and H39|l. Thus, 



(Bl) 



(B2) 



(B3) 



The hopping amplitude is somewhat reduced by the band narrowing factor. 
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